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We study macroscopic quantum tunneling of the magnetic moment in a single-domain particle 
placed above the surface of a superconductor. Such a setup allows one to manipulate the height of 
the energy barrier, preserving the degeneracy of the ground state. The tunneling amplitude and the 
effect of the dissipation in the superconductor are computed. 
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Tunneling of the magnetic moment in nanoparticles 
and molecular clusters has been intensively studied the- 
oretically and experimentally in the last decade jl| . The 
interest in this problem is two-fold. Firstly, magnetic 
tunneling reveals itself at a quasiclassical level, that is, 
in situations where all three components of the magnetic 
moment, M, can be rather accurately determined by a 
macroscopic measurement. The interaction of M with 
microscopic degrees of freedom makes this problem one 
of tunneling with dissipation Secondly, tunneling of 
the magnetic moment changes the magnetic properties of 
small magnets, with potential implications for the data- 
storage technology. It also adds nanomagnets to the list 
of candidates for qubits - the elements of quantum com- 
puters. 

In zero magnetic field the magnetic state of a classical 
magnet is degenerate with respect to — M due to 
time-reversal symmetry. In nanoparticles the | |> and 
| l> minima of the energy are separated by a barrier, 
U , due to the magnetic anisotropy. The thermal rate 
of switching between the two classical states is propor- 
tional to exp(—U/T). At high temperatures, when the 
thermal rate is high, the particle is in the superparamag- 
netic regime. At low temperature, as far as the thermal 
rate is concerned, the magnetic moment should freeze 
along one of the anisotropy directions. In particles of 
sufficiently small size, however, even at T = the mag- 
netic moment can switch due to quantum tunneling. The 
quantum switching rate, T, scales with the total spin S 
of the nanoparticle according |ri|-[l0| to In Toe — S. If the 
switching time, r , is small compared to the measuring 
time, the nanoparticle remains superparamagnetic in the 
limit of T = 0. If, in addition to that, the interaction of 
M with microscopic degrees of freedom (phonons, itiner- 
ant electrons, nuclear spins, etc.) is small, the nanopar- 
ticle, during a certain decoherence time, can exist in a 
coherent quantum superposition of two classical states, 
V&o = | T> +1 !>■ I n that state the probability that the 
moment of the particle has a certain orientation oscillates 



in time as cos(2ri). 

For a magnetic particle of considerable size to be in 
the quantum superparamagnetic (not necessarily coher- 
ent) regime, the energy barrier between the | j> and 
| !> states must be made sufficiently small. One way 
to decrease the barrier is to use an external magnetic 
field. This method has been used in experiments on mag- 
netic particles and molecular clusters performed to date 
It has a clear drawback if one attempts to cre- 
ate a coherent superposition of the | j> and | J.> states. 
Namely, the external field, unless it is applied exactly 
parpendicular to the anisotropy axis, removes the degen- 
eracy between the | f > and | |> states. In this Letter, 
we suggest a method of controlling the barrier without 
breaking the degeneracy. This method is illustrated in 
Fig. 1. The nanoparticle is placed above the surface 
of a superconductor at a variable distance controlled by, 
e.g., a piezoelectric layer or holder. The current induced 
in the superconductor creates the magnetic image of the 
nanoparticle. The interaction between the nanoparticle 
and the superconductor is then equivalent to the dipolc 
interaction between the nanoparticle and its image. The 
reduction of the barrier is similar to the one from the 
external magnetic field applied opposite to M. How- 
ever, contrary to the situation with the external field, 
the system (the nanoparticle plus the superconductor) 
is now degenerate with respect to M. It should 

be emphasized that such a degeneracy is a very general 
property of the system that is independent of the shape 
of the particle and the landscape of the superconducting 
surface. It is rooted in the time-reversal symmetry of the 
system in the absence of the field. The tunneling rate for 
the situation shown in Fig. 1 will be computed below. 

A high tunneling rate does not automatically provide 
the coherent superposition of quasiclassical states. Dif- 
ferent mechanisms of decoherence due to interactions in- 
side the nanoparticle have been worked out in recent 
years pd|-p^|. At low temperature, in the absence of 
nuclear spins and itinerant electrons, the effect of dissi- 
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pation on the tunneling rate can be very small |]. De- 
coherence is a more subtle issue. Generally speaking, 
macroscopic quantum coherence (MQC), that is cos(2rt) 
oscillations of M, occur only if the decohering interac- 
tions are small compared to the tunnel splitting of the 
ground state, HT. We will show that this condition can 
be satisfied at least as far as the interaction between the 
nanoparticle and the superconductor is concerned. This 
question is of interest also in a more general context of 
measuring the rotation of a mesoscopic spin with the help 
of a superconducting device. Indeed, most of the exper- 
iments on individual nanoparticles used SQUIDs. Al- 
though the problem studied here is different from those 
experimental situations, some of the ideas should also 
apply in those cases. 

Let us consider tunneling in a nanoparticle above the 
fiat surface of a superconductor PH| , as shown in Fig. 1, 
in the absence of dissipation. The external magnetic field 
will be considered zero throughout this paper. To make 
the classical electrodynamics of the problem less cum- 
bersome, we will make certain simplifying assumptions 
about the superconductor, the shape of the particle and 
its magnetic anisotropy. None of them is important and 
any generalization can be studied along the same lines. 
First, we shall assume that the superconductor is in the 
Meissner regime; that is, the magnetic field at the surface 
of the superconductor does not exceed H c i . We shall also 
assume that the characteristic geometrical dimensions of 
the problem, the size of the particle and its distance to 
the surface, are large compared to the penetration depth 
Xl. In that case, it is a good approximation to take the 
field of the particle at the surface of the superconductor 
to be parallel to that surface. The effect of the supercon- 
ductor on the particle is then equivalent to the magnetic 
dipole interaction with the image shown in Fig. 1. 

Next we assume that the particle is of ellipsoidal shape 
(that is, uniformly magnetized) with crystal fields either 
small or dominated by the single-ion anisotropy. The to- 
tal energy of the magnetic anisotropy of such a particle 
must be quadratic in the magnetization |15|] , 
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E an = yN ik MiM k 



(1) 



where M is the total magnetic moment of the particle, 
V is its volume, and the tensor Nik includes both the 
demagnetizing effect (that is, shape anisotropy) and the 
magnetocrystalline anisotropy. In a ferromagnetic par- 
ticle, M is proportional to V, while is independent 
of V. The factor I/ -1 in Eq. (1) is, therefore, needed 
to provide the correct linear scaling of E an with V. We 
shall assume that the principal axes of coincide with 
the coordinate axes in Fig. 1; with Z being the easy 
magnetization direction and N xx > N yy > N zz . The 
magnitude of M is assumed, as usual, to be formed by 
a strong exchange interaction and, thus, independent of 
the orientation. That is, Ml + M 2 + M 2 = M 2 = const. 



The energy of the magnetic anisotropy then becomes 
1 



E an = -{(3 x M' x -(3 z Mi) 



(2) 



where [3 X , (3 Z are positive dimensionless coefficients of or- 
der unity. 

Eq. (2) describes a magnet having a Y Z easy magne- 
tization plane with Z being the easy axis in that plane. 
The two degenerate minima of Eq. (2) correspond to M 
looking along and opposite to the Z axis. In the Meiss- 
ner state of the superconductor, the magnetic field of the 
particle induces superconducting currents whose field is 
equivalent to the field of the image shown in Fig. 1. 
As the particle moves closer to the superconductor, the 
interaction between the particle and its image increases 
and the barrier between the two equilibrium orientations 
of M decreases. The magnetic moment of the particle, 
M, and the moment of the image, m, are related through 



M x = m. 



My = m y , M z = 



(3) 



With reasonable accuracy the energy of the magnetic 
dipole interaction between the particle and its image is 
given by 



Ei n t 



Aim - 3(n-M)(n-m)] 

W 3 



(4) 



where d is the distance from the (center of the) particle 
to the surface of the superconductor. With the help of 
relations (4), one obtains |l4[ ] 



Ei n t — 



Ml 
(2d)* 



The total energy of the system, E = E a 
becomes 

E=^(-p z eMl+(3 x M 2 x ) , 



(5) 



E int then 



(6) 



where we have introduced e — 1 — V/ f3 z (2d) 3 . 

According to Eq. (6) the energy barrier between the 
degenerate | f > and | |> states of the particle is given by 
U = (3 z eM 2 V, where M a ~ M/V is the volume magneti- 
zation of the particle. Our main idea is to manipulate d 
in such a way that e and, consequently, U become small 
enough to provide a significant tunneling rate. Note that 
most deviations from the simplifying assumptions made 
above will renormalize f3 x . z and d in Eq. (6) but will 
not change the form of the total energy. It is convenient 
to introduce the total dimensionless spin of the particle, 
S = M./hj (7 being the gyromagnetic ratio) and two 
characteristic frequencies: 



u\\ = z eyM o , lj± = (3 x jM a 
The total energy can be then written as 



(7) 
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E= ~[-hu\\g + huxS%] 



(8) 



The corresponding tunneling problem has been studied 
by a number of authors [[t|-^0[. In the limit of W|| << luj_, 
i.e. at e << 1, the tunneling rate at T = in the absence 
of dissipation is given by T = A exp(—B ) with 
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(9) 



So far, we have neglected non-dissipative terms in the 
total energy that come from the superconductor. One 
such term is the kinetic energy of the Cooper pairs, 



E, 



n„mvt 



(10) 



where n s is the concentration of superconducting elec- 
trons, m is their mass, and v s = j s /en s is their drift 
velocity expressed in terms of the superconducting cur- 
rent j s - Eq. (10) can be written as 



_ _ 4ir\ 2 L 

E sc l — 5 — 



d 3 rf s 



(11) 



where = mc 2 /Aire 2 n s is the London penetration 
depth. The superconducting current is concentrated near 
the surface, resulting in the surface current [lq] 
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nxB(r) 



(12) 



where n is the unit vector in the Z direction and B(r) is 
the sum of the dipole fields of the magnetic particle and 
its image at z = 0. A somewhat tedious but straightfor- 
ward calculation then gives 



E scl = ^ [ d 2 rB 2 (r) 

87T J 



const 



3A L 
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Mi 



(13) 



The contribution of the kinetic energy of the supercon- 
ducting electrons to the total energy is small if d » X^. 
Even at c^Al, however, it reduces to the renormaliza- 
tion of the d dependence of the interaction between the 
magnetic particle and the superconductor. Since we are 
interested in d close to the critical value at which the 
barrier becomes zero, the form of the Hamiltonian, Eq. 
(8), and the expressions, Eq. (9), for the tunneling rate 
remain unaffected by that renormalization. 

The next non-dissipative term in the energy is the in- 
ertia coming from the energy of the electric field, 



d A r 
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(14) 



where A = — (47rA|/c)j s is the vector potential in the su- 
perconductor. Then, similarly to the previous case, one 
obtains 



E SC 2 — 



167TC 2 



3A ^ (M 2 



Ml) 



(15) 



To estimate the effect of this inertia term on tunneling, 
it should be compared, at the instanton frequency |l| 
Lj inst = (cjmwj,) 1 / 2 , with E of Eq. (8). The ratio of the 
energies is 



E sc2 /E~(\ L /d)\l/d){l/\ lnst f 



(16) 



where we have introduced / = V 1 ^ and \i ns t — c/^inst- 
Even if the experimental values of Z,e? and Al do not dif- 
fer in order of magnitude, Xi ns t can hardly be less than 
1 cm, making the ratio of energies in Eq. (16) negligible 
for nanoparticles used in tunneling experiments. We may 
then conclude that Eq. (9) gives a good estimate of the 
tunneling rate in the absence of dissipation. 

Spin tunneling with dissipation due to the interac- 
tion of M with microscopic degrees of freedom inside the 
nanoparticle has been intensively studied Inter- 
actions with phonons, magnons, nuclear spins, etc. have 
been considered. Here we will study the mechanism of 
dissipation specific to our problem: the interaction of 
the magnetic moment with normal quasiparticles in the 
superconductor. This analysis may also be relevant to 
experiments on spin tunneling performed using SQUIDs. 

We begin with the derivation of the energy dissipation 
in the superconductor due to the rotation of M. With 



the help of the relations E = 
one obtains 



l dK/dt and j = ctE, 
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d 3 ra(t)A 2 



(17) 



where a is the conductivity due to quasiparticles. 
Its Fourier transform can be approximated as er w = 
e 2 n q /m(v + iui), where n q is the quasiparticle concen- 
tration and v is their scattering rate. The latter is typ- 
ically 2-3 orders of magnitude greater than the instan- 
ton frequency for the tunneling of M. Consequently, 
the time dependence of a in Eq. (17) can be neglected 
and one can use a = e 2 n q (T) jmv. Equation (17) 
then becomes similar to Eq. (14) and the same argu- 
ment gives the surface integral that is proportional to 
(M 2 + M 2 )oz [9 2 + (9 2 + 2 )sin 2 6>]. Due to the smallness 
of e, the hard-axis anisotropy in Eq. (8) is small com- 
pared to the easy-axis anisotropy. Under this condition, 
quasiclassical trajectories of M must be close to the easy 
plane. This means 0«7r/2, while for the tunneling tra- 
jectory changes from to tt. A more rigorous analysis 
shows that <j) 2 ~e8 2 . Thus, with good accuracy, 



Q 



oA| 
2c 2 
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3A L M 2 
16vd 4 
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(18) 
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If Eq. (18) were quadratic in 9, it could be inter- 
preted as linear dissipation with a friction coefficient 
r\ = 3\LM 2 n q /16vd 4 n s . This would allow the Caldeira- 
Leggett approach to tunneling with dissipation. The 
dissipation in the rotation of M due to its interaction 
with quasiparticlcs is nonlinear in 9, however. Neverthe- 
less, since we only want to obtain an estimate of the effect 
of dissipation on tunneling, and because the sin 2 # term 
in Eq. (18) can hardly change this effect significantly, 
we shall go ahead and estimate the effective Caldcira- 
Leggett action as 



I C L = 



4tt 



dr' 



(It 



[9(r)-9(T')f 
(r- 



r A2 



(19) 



Note that the double-integral in Eq. (19) is dimension- 
less. 

The measure of the effect of dissipation on tunneling 
is the ratio Icl/^B where T%B is the effective action 
in the absence of dissipation, with B given by Eq. (9). 
After simple algebra, one obtains 



ICL 37T / ) 

Wl28</e \; 



d 



1 M 



(20) 



In a typical experiment, one should expect Xl < d and 
I < d, while the ratio ^Mq/v may hardly exceed 10~ 2 . 
Consequently, even at T^T C , when n q ^n s , the effect of 
the dissipation on the tunneling rate may become visible 
only at e less than 10 -5 . 

A more rigorous approach along the lines of Ref. flrj ] 
shows that the ratio n q /n s in equations (19)-(21) should 
be replaced by the ratio of the coherence factors. Similar 
to other dissipation problems due to quasiparticles [fl6|| , 
its effect on tunneling must have a maximum at T slightly 
lower than T c . Even at the maximum this effect should 
still be small. Coherence is a more subtle issue. It may 
be destroyed by a dissipative environment even if the ef- 
fect of dissipation on tunneling is weak. The mechanism 
of dissipation discussed above goes down with temper- 
ature as exp(— A/T), where A is the superconducting 
gap. Consequently, preserving coherent oscillations of M 
between | |> and | !> states requires T « A. 
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FIG. 1. Magnetic particle above the surface of a supercon- 
ductor. The effect of the superconductor on the particle is 
equivalent to the field from the image shown in the figure. 
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